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Abstract 

We study the differential branching ratio, branching ratio and the forward-backward asym- 
metry for the exclusive B s — > 4> ^ + ^~decay in the two Higgs doublet model. We analyze the 
dependencies of these quantities on the model parameters and show that these observables are 
highly sensitive to new physics and hence may provide powerful probe of the SM and beyond. 

1 Introduction 

The analysis of flavor-changing neutral current (FCNC) decays is one of the most promising re- 
search areas in particle physics from both theoretical and experimental sides. The rare B-meson 
decays induced by FCNC b — ► s transition, have received a special attention since their investi- 
gation opens up the possibility of a more precise determination of fundamental parameters of the 
standard model (SM), such as the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, 
the leptonic decay constants etc.. In addition, the studies on the rare B-meson decays open a 
window to investigate the physics beyond the SM, such as the two Higgs doublet model (2HDM), 
Minimal Supersymmetric extension of the SM (MSSM)[[lj], etc. to test these models and make 
estimates about their free parameters. 

The difficulties present in the experimental investigation of the inclusive decays stimulate the 
study of the exclusive decays. There exist now upper limits on the branching ratios of B° — > 
K°*n + fi- and B + -> K+n+fj,', given by CDF collaboration [g] 

BR(B° K°*/i + n-) < 4.0 x 1(T 6 
BR(B + K + fi + n-) < 5.2 x 1(T 6 . 

With these measured upper limits and also the recent measurement of the branching ratio of B — > 
K£+£~ with £ = e,n, 

BR(B -> Kt*T) = (0.75±g;|f ± 0.09) x 10~ 6 , 
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atKEK@|,the processes B — > (if, ET*)^ + £ have received great interest so that their theoretical 
calculation has been the subject of many investigations in the SM and beyond [Q]-[ 1 8 ]. Along 
this line, the exclusive decays induced by b — > s£ + £~ transition, like B s — ► <f)l + £~ , become also 
attractive since the SM predicts a relatively larger branching ratio, which is measurable in the near 
future. On the other hand, the theoretical analysis of exclusive decays is more complicated and 
contains substantial uncertainties due to the hadronic form factors, which need non-perturbative 
methods for their calculations. One of the methods that can be used to calculate the hadronic 
matrix elements is the three parameter fit of the light-cone QCD sum rule. In ref . Jl9|] , the form 
factors for b — > s£ + £~ induced exclusive B s — ► 4>£ + £~ decay have been calculated in the frame- 
work of this method. More recently, this process has been also investigated in a different context, 
namely using a so-called Constituent Quark-Meson model (CQM) based on quark-meson inter- 
actions [2C]. 

In this paper we investigate the exclusive B s — ► £ + £~ decay in the framework of the 2HDM. 
The transitions where B-meson decays into a vector meson, like K* fl]-JT8|] and p [21 ]-[^], have 
been extensively studied in the literature both in the SM and beyond. We calculate the dependence 
of the branching ratio (BR) and the forward backward asymmetry (Aps) of the exclusive B s — > 
4> ^ + ^~decay on the model parameters in the model I, II and III versions of the 2HDM and show 
that for some values of these parameters especially model III gives significance contributions to 
the BR and A FB . 

The paper is organized as follows: In section ^J, after we present the theoretical framework 
of the 2HDMs and the leading order QCD corrected effective Hamiltonian for the process b — ► 
s£ + £~ , we calculate the differential BR and the Afb of the exclusive B s — > <p ^ + ^~decay. The 
H section is devoted to the numerical analysis and the discussions. 



2 The theoretical framework 

Before presenting the details of our calculations, we would like to summarize the main essential 
points of the 2HDM, which is one of the most popular extensions of the SM. The 2HDM has two 
complex Higgs doublets instead of only one in the SM. In general, the 2HDM possesses tree- 



level FCNCs that can be avoided by imposing an ad hoc discrete symmetry [27]. As a result, 
there appear two different choices, namely model I and II, depending on whether up-type and 
down-type quarks couple to the same or two different Higgs doublets, respectively. Model II has 
been more attractive since its Higgs sector is the same as the Higgs sector in the supersymmetric 
models. Physical content of the Higgs sector includes three neutral Higgs bosons H°, h° and 
^4°, and a pair of charged Higgs bosons H^. In these models the interaction vertices of the 
Higgs bosons and fermions depend on the ratio tan (3 = v\/vz, where v\ and v% are the vacuum 
expectation values of the first and the second Higgs doublet respectively, and it is a free parameter 
in the model. The constraints on tan are usually obtained from B — B, K — K mixing, b — > 
decay width, semileptonic decay b — > ctv and is given by [ |28| ] 

0.7 < tan/? < 0.52(^1^), (1) 

and the lower bound m H ± > 200 GeV has also been given in 



In a more general 2HDM, namely model III [29, 30], no discrete symmetry is imposed and 
there appear FCNC naturally at the tree level. We note that in model III, FCNC receiving con- 
tributions from the first two generations are highly suppressed, which is confirmed by the low 
energy experiments. As for those involving the third generation, it is possible to impose some 
restrictions on them with the existing experimental results. Since the popular models I and II are 
special cases of the more general model III, we will give here the details of model III only. 
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The Yukawa Lagrangian in this general case is written as 



-D, 



£y = vYjQiL^iUjR + ri^QiL^DjR + ^'QiL^UjR + ^QiL^DjR + h.c. , (2) 



where i,j are family indices of quarks , L and R denote chiral projections L(R) = 1/2(1 =F 
75)> ipm form = 1,2, are the two scalar doublets, Qn are quark doublets, UjR, Djr are the 
corresponding quark singlets, r^- D and £^- D are the matrices of the Yukawa couplings. We can 
choose two scalar doublets ipi and if2 in the following form 



<Pi 



1 

72 




V2H+ 
Hi + iH 2 



(3) 



with the vacuum expectation values, 

< ipi >= 




< lf 2 > = 



(4) 



so that the first doublet ipi is the same as the one in the SM, while the second doublet contains all 
the new particles. Further, we take H\ and H2 as the mass eigenstates h° and ^4°, respectively. 

After the rotation that diagonalizes the quark mass eigenstates, the part of the Lagrangian that 
is responsible for the FCNC at the tree level looks like 



C y ,fc = -H^U[Vckm ^ v ckm L] V , 



(5) 



where U(T>) represents the mass eigenstates of up (down) type quarks. In this work, we adopt the 
following redefinition of the Yukawa couplings: 



(U,D 
SAT 



/ 4Gf £U,d 
V2 ^ N ' ij ' 



(6) 



Next step is to calculate the matrix elements for the inclusive b —> s£ + l~ decay in model III 
including the QCD corrections. For this, the effective Hamiltonian method provides a powerful 
framework. The procedure is to match the full theory with the effective theory, which is obtained 
by integrating out the heavy degrees of freedom, i.e., t quark, W^, H^, h° and H° in our case, at 
high scale fi = mw, and then calculate the Wilson coefficients at the lower scale \i ~ C(m&) us- 
ing the renormalization group equations. Following these steps above, one can obtain the effective 
Hamiltonian governing the b — > s£ + l~ transitions, in model III in terms of a set of operators 



n 



eff 



4G F 

71" 



10 



VtbV t *A aMOifa) 



(V) 



Here, 0\ and O2 are the current-current operators, Oz,...,0§ are usually named as the QCD pen- 
guin operators, O7 and 0% are the magnetic penguin operators and Og and 0\q are the semilep- 
tonic electroweak penguin operators. Cj(/i) are Wilson coefficients renormalized at the scale p,. 



The operator basis in the 2HDM for our process can be found in pl| , p2[ ]. 

Denoting the Wilson coefficients for the relevant process in the SM with Cf M (m w ) and 
the additional charged Higgs contributions with Cf(mw) » we have the initial values given by 



[31, 33] 



(mw) 




1 



3 



(jSM 



(m w ) 



3x? 



Cl M {m w ) 



r>SM 
u 9 



°10 



(m w ) 



4(s t - l) 4 
3s 2 
4(x* - l) 4 
1 



2x 2 , -8s 3 
* lnx t + 



5xf + 7xt 



24(x t - l) 3 
J + 5s? 

lnx t + 



s? + 5xf + 2xt 



8(x t - If 



. 1 - 4sin 2 6»vy , _.. , 4 
— 37-5 (x t + . 2 - C(x t ) - D(x t ) + -, , 



1 



sin" Uw 



-(B(x t )-C(x t )) , 



and 



(m w ) 
C?{m w ) 
C£{m w ) 

Ci {m w ) 



0, 

Y 2 F 1 (y t ) + XYF 2 {y t ) , 
Y 2 G 1 {y t ) + XYG 2 {yt) 
Y 2 H 1 {y t ) , 
Y 2 L x {y t ) , 



(8) 



where 



x t 
X 
Y 



m 



Vt 



mi 



m 



1 

m b 
1 

m t 



FD i FD *td 
$N,bb ?7V,db"^ 

37 i 37 



(9) 



Note that the results for model I and II can be obtained from model III by the following substitu- 
tions: 



Y 
Y 



cot p , XY 

cot (3 , XY — * 1 for model II . 
The explicit forms of the functions -F\(2) {yt), ^1(2) (yt)> Hi (yt) and Li (y t ) in Eq.(|8|) are given as 

yt(7-5yt w - 2 ' " 2 ' 



cot /? /or model I 



Fi{y t ) 

Gi(y t ) 
G 2 (yt) 
Hi(y t ) 



72(y t - 1)3 

yt(% 



8yf) y (3y t - 2) 

+ — — In yt 



12(yt - 1) 4 



3) y f (-3y t + 2) 

+ — rr^ Inyt 



12(^-1)2 ' 6(y t -lf 
yt(-y 2 + 5y t + 2) | -y 2 

24(y t - l) 3 
2/t(y* - 3) , yt 



Kyt - 1)'- 



ln yt , 



+ 



4(y t - l) 2 2(y t 
1 — 4sm 2 #jy xy t 
sin 2 6\Y 8 



1)3 
1 



In yt , 



1 



Vt 



yt 



Li(yt) 



A7y 2 - 79y t + 38 
108(y t - l) 3 



xyt 



- 1 (yt - i) 2 

3yf - 6y t + 4 



In yt 



sin 2 6\y 8 



1 



yt - i (y* - 1) 



i) 

2 ln J/t 



(10) 
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Finally, the initial values of the coefficients in the model III are 

Cf HDM (m w ) = C! M (m w ) + C*(m w ). (11) 

Using these initial values, we can calculate the coefficients Cf H (fj) at any lower scale in the 
effective theory with five quarks, namely u, c, d, s, b similar to the SM case J3T|]-[|34]]. 

The Wilson coefficients playing the essential role in this process are Cj HDM '(//), Cg HDM (//), 
Cfp D (fi). For completeness, we also give their explicit expressions 

c^(») = cr DM {») + Q d (ci HDM (n) + n c cl HDM M) , 

where the leading order QCD corrected Wilson coefficient C 7 ° ,2HDM {[i) is given by [ ]34[ [3l| , 32] 

c LO,2HDM {fi) = ^16/23 C72 HDM (mw) + (8/ 3 )( ^14/23_^16/23 )c 2HDM (mw) 

+ Ci HDM (m w )J2hiV ai , (12) 

i=l 

and r] = a s {mw)l 'a s (/x), hi and a, are the numbers which appear during the evaluation [34]. 

Cg (fi) contains a perturbative part and a part coming from LD effects due to conversion of 
the real cc into lepton pair £ + £~: 

C 9 e// (/x) = <T rt ( M ) + Y reson {s) , (13) 

where 

+ h(z, s)[3Ci(/x) + C 2 (At) + 3C 3 (/*) + C 4 (m) + 3C 5 (^) + C 6 (m)] 

- s) (4C 3 (m) + 4C 4 (ax) + 3C b (m) + C 6 (/i)) 

- i/ i (0, S )[C 3 (/x) + 3C 4 ( / u)] (14) 

+ |(3C 3 (/i) + C4(M) + 3C 5 (/i) + C 6 (M)) , 

and z = m c /m b . The functions s) arises from the one loop contributions of the four quark 
operators Oi,...,Oq and their explicit forms can be found in [[32|, Q 

It is possible to parametrize the resonance cc contribution Y reson (s) in Eq.([TJ) using a Breit- 
Wigner shape with normalizations fixed by data which is given by [ ^ ] 

Y (si - 3 k V ^(^^'W, 

a em v ^p % Q ~ m-Vi + im Vi T Vi 

x [(3C7 1 (//) + C 2 (M) + 3C , 3(M) + C 4 (M) + 3C7B(At)+C , 6(M))]- (15) 

The phenomenological parameter k in Eq. ( p3| ) is taken as 2.3 so as to reproduce the correct value 
of the branching ratio BR(B -» J/tp X -» XtJt) = BR(B -» J/^ X)BR{J/^ -> X££). 

Neglecting the mass of the s quark, the effective short distance Hamiltonian for the 5 — > s^ + ^ _ 
decay leads to the QCD corrected matrix element: 

M = ^^V tb V t * s \cl ff {m b ) s 7m (1 - lb )bhH + C w (m b ) s^(l - j 5 )b£^ l5 £ 
- 2C e 7 ff (m b ) ^Ma^il + ^bl^, 

(16) 

where q is the momentum transfer. 
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2.1 The exclusive 



£+£ decay in the 2HDM 



In this section we calculate the BR and the Ap B of the B s — ► <j> £ + £ decay. In order to find 
these physically measurable quantities at hadronic level, we need the following matrix elements 

< 0(p0,e)|s7 M (l - m)b\B(p B ) > and < 4>(p^, e)\si<T IMV q v (l + j 5 )b\B(p B ) >, which can be 



parametrized in terms of form factors. Using the parametrization of the form factors as in [ 12], 
we get the following expression for the matrix element of the B s — ► (f> l + 1~ decay: 



M L 



Gfci 



i^VtbVtlU^Ae^e^pye + iBe* - iC (p B + p^s* q) - iD(e*q) qil } 



+ h^d[2Ee^ Xa e* u P y B + iFe^ - iG{e*q)(p B - iH{e*q)q^\ 



(17) 



where 



.4 

B 

C 

D 

E 
F 
G 

H 



a 



eff 



V 



m B + m. 



(m B + m*) (cl ff A l + ^L{ mB - m )C 7 e// T^ 



n e ff 



An 



m B + nicf, q 



m B -m- 



-T 3 , 



2C^^(A 3 -A )-4C^^T 3 , 
q z q z 

r v 

WO ; , 

m B + m<f, 

Ci (m B + m^Ai, 
A 2 



C 



10 



2C W ^(A 3 -A ), 



(18) 



Here Aq, A\, A 2 , A 3 , V, T\, T 2 and T 3 are the relevant form factors. 



The matrix element in Eq.(|17|) leads to the following differential decay rate Q17Q 

dT _ a 2 G 2 F m B 2 



ds 



2 12 TT 5 



WtbV ts \ 



'Xs v A„ 



where 

Aw, = 



+ 



-A^m|s((3 - v 2 )\A\ 2 + 2v 2 \E\ 2 ) + — X^m% 



1 

>'<P 



(19) 



-X^m B (3 - v 2 )\C\ 2 + m B s z (l - v z )\H\ 



+ o t( 3 - y2 ) s (! - v 2 )]Re[F G*] - 2 s (1 - v 2 )Re[F H*} 



+ 2m|s(l - t>)(1 - v 2 )Re[G H*] + -(3 - v 2 )W<pRe[B C*] 

3 



3r, 



-m B 



(A + 12r s)(3 - ?; 2 )|5| 2 + A m|[A (3 - v 2 ) - 3s(s - 2 r<t> - 2)(1 - w 2 )]|Gp 



+ (A (/) (3-7; 2 ) + 24r^t; 2 )| J F| : 



(20) 
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QCDSR 


CQM 






(Jj p 


up 




Li 


up 




0.30 


0.87 


-0.06 


0.59 


-0.11 


0.49 


At* 


0.26 


1.55 


0.51 


0.73 


0.78 


-0.52 


yB-** 


0.43 


1.75 


0.74 


0.20 


0.65 


0.96 


rpB-><l> 


0.35 


1.82 


0.83 


0.21 


0.78 


0.07 


rpB-Hf) 


0.35 


0.70 


-0.32 


0.21 


0.85 


12.9 


rpB-*4> 

J 3 


0.26 


1.52 


0.38 


0.18 


0.62 


-0.88 



Table 1 : The values of the parameters in Eq. ( |2q ) for the various form factors of the transition B — > </> 
calculated in the light-cone QCD sum rule approach (QCDSR) [[T9p and using a Constituent Quark- 
Meson model (CQM) [JZtJ. 



where s = q 2 /m 2 B , = m^/m|, v = t = mi/m B , = + (s - l) 2 - 2r 4> {s + 1) 

and W(f, = — 1 + + s. Here, z = cos 6, where 8 is the angle between the three-momentum of 
the l~ lepton and that of the B s -meson in the center of mass frame of the dileptons l + t~ . 

The Aps is another observable that can give more precise information at hadronic level. We 
write its definition as given by 

Afb{s) = Y d /r n (2D 

where T is the total decay rate. The Aps for the B s — ► ^ + £~decay is calculated to be 

A FB = /" ds 8mBA<^u 2 s(i?e[S J5*] + J2e[>l F*]) J J dsyfx^ v A . (22) 



3 Numerical results and discussion 



In this section we present the numerical analysis of the exclusive B s 
2HDMs. The input parameters we used in this analysis are as follows: 



"decay in the 



m Bs = 5.28 GeV , m b = 4.8 GeV , m c = 1.4 GeV , m M = 0.105 GeV , 
m T = 1.77 GeV, = 1.02 GeV ,m H ± = 400 GeV , \V tb V t *\ = 0.04, 
a' 1 = 129 , G F = 1.17 x 10~ 5 Gel^ 2 , r Sa 



1.54 x 10~ 12 s. 



(23) 



The masses of the charged Higgs, m H ±, the Yukawa couplings (£^' D ) and the ratio of the 
vacuum expectation values of the two Higgs doublets, tan j3, remain as free parameters of the 
model. The restrictions on m H ±, and tan j3 have been akeady discussed in section 0. For Yukawa 



couplings , we use the restrictions coming from CLEO data Q36|], 



BR{B X s 7 ) 



(3.15 ±0.35 ±0.32)10" 



(24) 



B° 



NAb 



B° mixing [33], p parameter []30t], and neutron electric-dipole moment [37], that yields 
' L> 1 0, where the indices i, j denote d and s quarks, and ^Ntc ^Ntt- 



and i L N ij 



Therefore, we take into account only the Yukawa couplings of b and t quarks, £ 



v 

N.tt 



> CaU- There 
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is also a restriction on the Wilson coefficient C^" from the BR of B 



X s 7 in Eq.(24) as follows 



[33], 



0.257 < \CZ ,f \ < 0.439. 



(25) 



In our numerical calculations for B s — ► <j> £ + £ decay, we use three parameter fit of the 



light-cone QCD sum rule [ |19| ] which can be written in the following form 



F(q 2 ) 



F(0) 



I -dp q 2 /m? B + b F (q 2 /m 2 B ) 2 



(26) 



where the values of the parameters F(0), ap and bp are given in Table ([]]). This table also contains 



the values of the same parameters calculated in the CQM [2C], which we use to calculate the SM 
predictions for the BR in table ^ The form factors Aq and A3 in Eq. (26) can be found from the 
following parametrization, 



A 
A 3 



A, 



T 3 g 2 



2m^ 



m B - 
2m^ 



A*. 



(27) 



We note that there are five possible resonances in the cc system that can contribute to the 
decay under consideration and to calculate them, we need to divide the integration region for q 2 
into three parts for £ = e, fx so that we have Am 2 < q 2 < (m^ 1 — 0.02) 2 and (m^ 1 + 0.02) 2 < 
q 2 < {m^ 2 — 0.02) 2 and {m^ 2 + 0.02) 2 < q 2 < (mp — m^) 2 , while for £ = r it takes the form 
given by 4m 2 < q 2 < (m^ 2 — 0.02) 2 and (m^ 2 + 0.02) 2 < q 2 < (tub — m^) 2 .Here, m^, 1 and 
m^p 2 are the masses of the first and the second resonances, respectively. The SM predictions for 
the integrated branching ratios for £ = e, /j,, r are presented in Table ^| for distinct q 2 regions and 
also the total contributions. 

In the following, we give the results of our calculations about the dependencies of the dif- 
ferential branching ratio (dBR/dq 2 ) and Apsiq 2 ) of the B s —> <fi ^ decay on the four 
momentum transfer q 2 , and also dependencies of the BR and App on the model parameters . The 
results are presented by a series of graphs, which are plotted for £ = \i and for the case of the 



ratio \r tb \ 



< 1. We do not present the results for r tb > 1 case, since BR for r t & > 1 



case indicates one-to-two orders of magnitude enhancement with respect to the SM results, which 



seems to conflict with the experimental data given by BELLE collaboration [ J39| ] for the inclusive 
B -> X s £ + £- decays : 



BR(B -> X s 



1 



(7.1 ± 1.6±H) x 10" 



(28) 



In Fig.([I]), we plot the dependence of the BR of the B s — > cj) ^ + /u~decay on tan j3 by taking 
mH± = AQQGeV . Here the solid (dashed) curve represents the Model I (II) prediction for the BR 
and the small dashed straight line is for the SM result. We see that BR decreases with increasing 
values of tan (5. For 1 < tan (5 < 4, it is possible to enhance the BR 15(13)% — 2(1.5)% in model 
I(II) compared to its value in the SM. However the larger values of tan j3 have been favoured 
by the recent experimental results [38]. We can therefore conclude that charged Higgs boson 
contributions calculated in the context of the model I and II versions of the 2HDM are not very 
sizable, i.e., when tan /3 > 5 it almost coincides with the SM result in model I, while in model II, 
its value remains slightly below the SM one. 
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Table 2: The SM predictions for the integrated branching ratios for £ = e, fi, r of the B s — > 
£ + £~decay for distinct g 2 regions and also the total contributions. Here, QCDSR (CQM) stands 
for the results calculated with the form factors of ref. [ |19| ] ([|20|]). 

We have the following prescription for the graphs we give in between Figs.(||)-(|6|): the regions 
between the solid curves represent the > case, the regions between the dashed curves 

represent the < case, and the small dashed curves are for the SM predictions. 

The dependence of the dBR/dq 2 on q 2 is given in Fig.(^) in Model III, including the long 
distance contributions. We see that dBR/dq 2 almost coincides with the SM result for Cy < 0, 
while for Cy > cases it is considerably enhanced. 

The dependence of the BR on one of the free parameters of Model III, |^ bb /m,b is shown in 
Fig.@. It is seen from this figure that BR is not very much sensitive to fjy bb/ m b> especially for 
large values of this parameter. We note that the SM and model III average results for BR when 
Cj ff > (Cf f < 0) are 1.91 x lO" 6 and 4.00 x lO" 6 (1.91 x 10~ 6 ). 

The dependence of the BR on the charged Higgs mass m H ± is presented in Fig.(||). The BR 
decreases as m H ± increases, except the C% * < case, which is insensitive to m H ± and very 
close to the SM prediction in magnitude. 

We present the dependence of the Afb{<1 2 ) on q 2 for the B s — > c/> /i + /i-decay in Fig.(Q). 
Here, the Afb(q 2 ) is enhanced for > while it almost coincides with the SM prediction 
for Cj ff > 0. 

In Fig.(Q), the dependence of the ipg on ^Nbb/ m b i s represented. We observe that Afb is 
not very sensitive to fjy bb/ m b > especially for large values of this parameter. The SM and model 

III average results for A FB when cf f > (Cf f < 0) are -0.23 and -0.31 (-0.21). 

In conclusion, we have investigated the physical observables, BR and Afb related to the 
exclusive B s — ► <j> £ + £-decay in the model I, II and III versions of the 2HDM. We have found 
that these observables are highly sensitive to new physics and hence provide powerful probe of 
the SM. 
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Figure 1: The dependence of the BR of the B s — > <p At + /i~decay on tan (3. Here the solid (dashed) 
curve represents the Model I (II) contribution while the small dashed straight line is for the SM case. 




Figure 2 : The dependence of the dBR/ dq 2 onq 2 . Here the region between the solid curves represents 
the dBR/dq 2 for C e 7 ff > 0, while the one between the dashed curves is for Cj ff < 0. The SM 
prediction is represented by the small dashed curve. Here we take £jy bh = 40 m b . 
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Figure 3: The dependence of the BR on ^ bb /m b . Here BR is restricted in the region between solid 
(dashed) curves for Cy > (C^ < 0). Small dashed straight line represents the SM prediction. 
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Figure 4: The dependence of the BR on m H ± for £jy bb = 40 m b . Here BR is restricted in the region 
between solid (dashed) curves for Cy > (Cy < 0). Small dashed straight line represents the SM 
prediction. 
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Figure 5: The dependence of the A FB (q 2 ) on q 2 . Here the region between the solid curves represents 
the A FB (q 2 ) for C e 7 ff > , while the one between the dashed curves is for the A FB {q 2 ) for C 7 e// < 0. 
The SM prediction is represented by the small dashed curve. Here we take £® bb = 40 m b . 
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Figure 6: The dependence of the A FB on £^66- 



(dashed) curves for C% * > (Cy 77 < 0). Small dashed straight line represents the SM prediction. 



Here A FB is restricted in the region between solid 
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